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1. Introduction 

This is the first of four proceedings providing an update of our calculation of Bk using im- 
proved staggered fermions (HYP-smeared valence quarks on the MILC asqtad lattices). Here we 
present results using fits based on SU(2) staggered chiral perturbation theory (SChPT), the method 
that leads to our most accurate results. We focus on the progress made since last year's lattice pro- 
ceedings, Ref. [JXJ] . The other proceedings present, respectively, results using fits based on SU(3) 
SChPT [Bj, a study of some sources of error [H], and a more detailed look at results from the 
ultrafine ensemble [Q]. 

Table [j] shows the present status of our running. In the last year, we increased the number 
of measurements on the C4 and SI ensembles (by factors of 10 and > 2, respectively), and added 
two new ensembles: F2 and Ul. The F2 ensemble allows us to further check the sea-quark mass 
dependence, while the Ul ensemble provides a fourth lattice spacing. We have also written a long 
article, Ref. [S], in which we explain both the SU(2) and SU(3) SChPT calculations leading to our 
fitting forms, describe our fitting methods, and present our full error budget. The results presented 
in Ref. [j|] are intermediate between those presented at last year's lattice conference [|l|] and those 
presented here. In particular, Ref. [|5|] does not include results from the Ul ensemble. 



a (fm) 


ami /am s 


geometry 


ID 


ens x meas 


B K Qi = 2 GeV) 


status 


0.12 


0.03/0.05 


20 3 x 64 


CI 


564 x 1 


0.557(14) 


old 


0.12 


0.02/0.05 


20 3 x 64 


C2 


486 x 1 


0.569(16) 


old 


0.12 


0.01/0.05 


20 3 x 64 


C3 


671 x 9 


0.565(5) 


old 


0.12 


0.01/0.05 


28 3 x 64 


C3-2 


275 x 8 


0.570(5) 


old 


0.12 


0.007/0.05 


20 3 x 64 


C4 


651 x 10 


0.562(5) 


update 


0.12 


0.005/0.05 


24 3 x 64 


C5 


509 x 1 


0.554(11) 


old 


0.09 


0.0062/0.031 


28 3 x 96 


Fl 


995 x 1 


0.544(12) 


old 


0.09 


0.0031/0.031 


40 3 x 96 


F2 


678 x 1 


0.547(10) 


new 


0.06 


0.0036/0.018 


48 3 x 144 


SI 


744x2 


0.539(7) 


update 


0.045 


0.0028/0.014 


64 3 x 192 


Ul 


305 x 1 


0.527(11) 


new 



Table 1: MILC asqtad ensembles used to calculate B K . Z?a:(NDR,2 GeV) is obtained using the SU(2) 
analysis with 4X3Y-NNLO fits (discussed in the text). 

2. SU(2) SChPT Analysis 

Our analysis makes essential use of SChPT [@, 0] in order to remove artifacts associated with 
taste breaking. The application of SChPT to Bk was worked out for SU(3) SChPT in Ref. [j|]. 
The extensions to our mixed-action set-up, and to SU(2) ChPT, were presented in Ref. [||]. We 
summarize the findings of the latter work in this section. 

SU(2) ChPT treats the strange quark as heavy, and expands in m\jm\, as well as in the usual 
ratio m\j K 1 ^ (with A x ~ 1 GeV) [l(J. It has been argued that this is a more reliable way of 
extrapolating kaon properties to the physical quark masses than using SU(3) ChPT, in which the 
strange quark is treated as light. Whether or not this is true in general, it turns out that, for our 
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application, there is an additional important advantage of SU(2) ChPT. This is that, at next-to- 
leading-order (NLO), the SU(2) SChPT expression contains no low-energy coefficients (LECs) 
arising from discretization or truncation errors. 1 This is not the case in SU(3) ChPT, and, as a 
result, the latter gives rise to cumbersome fitting forms with many parameters. In SU(2) SChPT, by 
contrast, one has the same number of parameters at NLO as for a fermion discretization with chiral 
symmetry, such as domain-wall fermions. We find ^ 

/th = d Q F + d\ -jg + d 2 ^j + d 3 -^- (2. 1) 

where the chiral logs reside in the function 

F ° = i + ^^IW + ^i-W^)- 2 ^))} (2-2) 



32^ 

<W> = ^ 2 \e(x 1 )+e(x P )+4e(x w )+4e(x A ) + 6£{x T 



16 L- - - J ' (23) 

£(X) = X [log(X//ig R )] , l(X) = -d£(X)/dX (2.4) 

Here X B (L B ) is the squared mass of the valence (sea) pion with taste B, which we know from 
our simulations or those of the MILC collaboration. It is only through these masses that the taste- 
breaking enters at NLO. The coefficients d, have an unknown dependence on r s = m s / 'Aqcd, and, 
in addition, do is a function of a 1 and q s 2 . 



We have included an analytic NNLO term in Eq. |2.1| , with coefficient d 2 , and we do fits both 
without and with this term, labeled respectively as NLO and NNLO. The latter fits are not full 
NNLO fits, but rather are used to gauge the errors arising from truncating ChPT. 

3. Fitting and Results 

Our lattice kaon is composed of valence (anti)quarks with masses m x and m y , which are to be 
extrapolated to m p d ys and mf ys , respectively. On each lattice we use 10 valence masses: 

am x , am y = am s x n/10 with n = 1,2,3, 10, 

where am s is the strange sea-quark mass (which lies fairly close to the physical strange mass). 
For our standard SU(2) fitting we choose the lowest 4 values for am x and the highest 3 values for 
amy, calling the resulting fits "4X3Y". This means that we use only 12 out of the possible 55 mass 
combinations, but by doing so we maintain the SU(2) condition m x /m y <C 1. 

We do the fitting in two stages. First, for fixed m y , we fit the m x (or, equivalently, Xp) depen- 
dence to the form (2T). We call this the "X-fit", examples of which are shown in Fig. [jj This fit 



allows us to extrapolate m x — > m j hys and to remove taste-breaking discretization and perturbative 



truncation errors by setting Xj = Xb = X phys in Eq. (2.2). We can also partially correct for the use 



of an unphysical mi by setting L/ to its physical value. The result of this procedure is shown in the 
figures by the red points. 



1 Truncation errors arise because we match the four-fermion operator appearing in Br to continuum regularization 
using one-loop perturbation theory. 
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Figure 1: 4X-NLO fit of Bk versus Xp for the C3 (left) and SI (right) ensembles. The red point gives the 
physical result, as discussed in the text. Here, we fix am y = 0.05 for C3 and am y = 0.018 for SI. 



In the second stage, we extrapolate the corrected points from the 3 values of m y to m!j >hys . The 
dependence here is analytic, and we use both linear and quadratic forms in these "Y-fits". The 
Y-fits are straightforward, and we refer to Ref. [||] for examples. Here we focus on the X-fits. 



ID (meas) 


do 


d\ 




Z 2 /d-o.f 


C3 (671 x 9) 
SI (744 x2) 


0.5602(34) 
0.4808(50) 


0.035(16) 
0.143(30) 




0.83(54) 
0.06(14) 


SI (513) 
SI (744 x 2) 


0.469(15) 
0.4834(93) 


0.33(21) 
0.09(13) 


-0.72(81) 
0.27(51) 


0.002(18) 
0.075(94) 



Table 2: Parameters of X-fits shown in Figs, [TJand [| 

The parameters of the fits shown in Fig. [I] are given in the first two rows of Table |. 2 We use 
uncorrected fits and thus expect # 2 /d.o.f <SC 1 for a good fit. The C3 data are unchanged from 
last year (Refs. [jl|, ||]) while the SI results have significantly improved statistics. The curvature in 
these NLO fit functions is entirely due to the chiral logarithms, and is consistent with our data on 
all ensembles. The convergence of the chiral expansion can be gauged from the difference between 
the values of do in the Table (which is the LO result) and the results in the figures. The ratio of 
NLO to LO terms is < 10% for all points used in the fits. Such satisfactory convergence is seen on 
all other ensembles [j|]. 

The rather large % 2 in the C3 fit may be indicative of the need to include NNLO terms, and 
indeed the % 2 of our NNLO fit is much reduced. The resulting extrapolated-corrected value of Bk 
(the red point) is, however, very similar in both fits. 

In Fig. |2[ we compare 4X-NLO fits on the Fl and F2 lattices. The former results are from last 
year, while the latter are new. These two ensembles differ only in the value of light sea-quark mass, 

2 For a given ensemble with fixed ami, ^3 i s absorbed into dy. 
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Figure 2: Bk versus Xp on the Fl (left) and F2 (right) ensembles, showing a 4X3Y-NNLO fit. Notation as 
in Fig. [j] Here, we set am y — 0.031 for both Fl and F2. 
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Figure 3: Bk versus Xp on SI ensemble for 513 configurations (left) and for 744 configurations x 2 mea- 
surements (right). The plots are for am y = 0.018. 4X3Y-NNLO fits are shown. Notation as in Fig. [j]. 



and we see that results are very similar on the two ensembles. This can be seen quantitatively from 
the resulting values of Bk(2 GeV) given in Table |I| The lack of dependence on ami confirms the 
result found on the coarse ensembles C1-C5. 

In Fig. ||, we show the impact on X-fits of improving the statistics on the SI ensemble. We 
compare 4X3Y-NNLO fits between last year's statistics and our new results. The corresponding fit 
parameters are given in Table ||[ 

Increasing the statistics makes the visual evidence of curvature more convincing and improves 
the determinations of the LECs. Comparing the SI results from NLO fits and NNLO fits we see 
from Table 0) that do and d\ are consistent, although the latter is poorly determined in both fits. 
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Fortunately, it is da which gives the dominant contribution to the physical Bk, and this is well 
determined by both fits. 

4. Continuum Extrapolation 

Repeating the above procedure on all ensembles yields the results for Bk quoted in Table [j]. 
Here we have run these results to a common renormalization scale, /I = 2 GeV. The next steps are 
to extrapolate to the continuum limit, and to estimate all sources of error. 

We do the continuum extrapolation using the results from the C3, Fl, SI and Ul lattices. These 
lattices all have the same ratio of sea quark masses, nn/m s = 1/5 and all have m s m^ hys . We have 
seen above that Bk is almost independent of the sea-quark mass, so the lack of exact matching of 
the sea-quark masses between these ensembles has a very small effect (and can be corrected for). 

The expected approach to the continuum limit is somewhat complicated. The dominant errors 
remaining in Bk are due to taste-conserving discretization errors and errors from the truncation of 

matching factors. The discretization errors have the form a 2 a s {\/a) n , with n = 0,1,2, Note 

that n = is allowed because we do not Symanzik-improve our operators. The truncation errors, 
by contrast, are proportional to a s {\/a) m with m = 2, 3, 

Since we cannot hope to disentangle these various dependencies, we adopt a pragmatic ap- 
proach. We fit to the form 

B K (a) =B K (a = 0)+bia 2 + b 2 a 4 , (4.1) 

with or without the quadratic term. This takes care of the dominant a 2 term, and the inclusion of 
the &2 is an approximate way of allowing for dependence such as a 2 a s . Clearly this fit does not 
treat the truncation error correctly. Instead, we treat this as a systematic error, as discussed below. 

In Fig. [j], we show linear and quadratic fits to the a 2 dependence of Bk- The fits agree well, and 
we choose the linear fit for our central value and statistical error. We quote the difference between 
the result on the Ul lattice and the continuum extrapolated value as an estimate of systematic error 
in the continuum extrapolation (due to our not accounting for a 2 a s terms). 

5. Error Budget and Conclusion 



cause 


error (%) 


memo 


status 


statistics 
matching factor 
discretization 
fitting (1) 
fitting (2) 
ami extrap 
am s extrap 
finite volume 

n 


1.4 

4.4 

0.10 

0.92 

0.08 

0.48 

0.5 

0.85 

0.14 


4X3Y-NNLO fit 
Aflg } (Ul) 

diff. of (Ul) anda = 
X-fit (C3) 
Y-fit (C3) 

diff. of (C3) and linear extrap 
constant vs linear extrap 
diff. of 20 3 (C3) and 28 3 (C3-2) 
ri error propagation 


update 
update 
update 

i 
© 
i 
@] 
i 





Table 3: Error budget for B K obtained using SU(2) SChPT fitting. 
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Figure 4: Continuum extrapolation of Br-(NDR,jII = 2 GeV). Linear and quadratic fits to a 2 (in fm 2 x 100) 
are shown. The results are from 4X3Y-NNLO SU(2) SChPT fits. 



The error budget for the SU(2) analysis is presented in Table |3| Several of the errors are 
as in Ref. [Q], and we refer to that work for explanations and discussion of these errors. 3 The 
errors which have changed since Ref. ^ are those due to statistics (reduced from 1.7%), matching 
(reduced from 5.5%) and discretization (reduced from 1.8%). Our estimate of the latter error, 
explained in the previous section, may be an underestimate, but in any case is dominated by the 
matching error. 

The matching error is estimated as follows. We assume that the dominant missing term in the 
perturbative matching factors is of size 1 x (Z^ (jU = I /a), so that the error in Bk is 



H 2) «4 1] x 



a,(l/a) 



2 



(5.1) 



where is the result using one-loop matching. This error will be reduced, but not eliminated, 

(2) 

by the continuum extrapolation. To be conservative, we take as the error the size of AB K on our 
smallest lattice spacing. The reduction in this error compared to Ref. ^ is simply due to our having 
another, smaller, lattice spacing. 

Combining the errors in the error budget, we find 

^ (NDR, n = 2 GeV) = 0.5260 ± 0.0073 ± 0.0244 

B K = B K {RGl) =0.720 ±0.010 ±0.033. °'~ ! 



Our result is in agreement with those obtained using valence DWF on either MILC flllQ or DWF 



lattices [[R], 12]. Our total error of 5% is somewhat larger than the 4.1% and 3.3% attained in 
the other calculations, the difference being mainly due to our use of one-loop, rather than non- 
perturbative, matching. 



3 See also the companion proceedings ^ for additional discussion of errors due to finite- volume effects and the 
dependence on sea-quark masses. 
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